Organic superconductors: the need to go beyond effective ^-filled band models 
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A frustrated, effective ^-filled band Hubbard-Heisenberg model has been proposed to describe 
the strongly dimerized charge-transfer solid families k-(ET)2X and Z[Pd(dmit)2]2- In addition to 
unconventional superconductivity these materials also exhibit antiferromagnetism, candidate spin- 
liquid phases, and in the case of Z=EtMe3P, a so-called valence-bond solid phase. We show that 
neither superconductivity nor the valence-bond solid phase occurs within the Hubbard-Heisenberg 
model, indicating that the effective ^-filled band model is unsuitable for these materials. 

PACS numbers: 71.10.Fd, 71.10.Hf, 74.20.Mn, 74.70.Kn 



Low-dimensional organic charge transfer solids (CTS) 
are being intensively studied because of their many un- 
usual competing and coexisting electronic phases. The 
most studied among them are probably the k-(ET) 2 X 
and Z[Pd(dmit)2]2 families, which, depending on the an- 
ion X - or cation Z + exhibit unconventional superconduc- 
tivity (SC), Neel antifcrromagnetic (AFM) order, charge 
ordering (CO), candidate quantum spin liquid (QSL) be- 
havior, and valence-bond solid (VBS) order The ap- 
parent similarity between these with the cuprate super- 
conductors have been noted by many investigators. The 
semiconductor-SC transition in the CTS occurs under 
the application of pressure at constant carrier density, 
which suggests that the transition is driven by a small 
modification of a particular parameter of an appropri- 
ate Hamiltonian. The key questions then are, what is 
the minimal model, and which is the parameter whose 
changes give the competing phases. 

Experimental observations appear to give a simple an- 
swer to these questions. The 2:1 (1:2) stoichiometry of 
k-(ET) 2 X (Z[Pd(dmit) 2 ] 2 ), indicates 0.5 hole (electron) 
per cationic (anionic) molecule. The crystal structures 
consist of dimers of molecules arranged in an anisotropic 
triangular lattice. The strong role of repulsive electron- 
electron (e-e) Coulomb interactions is indicated by the 
observation of Neel AFM in k-(ET) 2 Cu[N(CN) 2 ]C1, a 
candidate QSL phase in k-(ET)9Cu 2 (CN)3, and a VBS 
phase in EtMe 3 P[Pd(dmit) 2 ] 2 [l(. Taken together, the 
effective ^-filled band Hubbard model, or some variant 
of it, now appears to be the appropriate minimal model, 
with the dimer units as the sites. The anisotropy of the 
triangular lattice, i.e., the degree of frustration, is the 
variable parameter that changes under external pressure 
or internal pressure effect caused by counter-ions with 
large size. Such a picture readily explains the observed 
commensurate AFM at large anisotropy. Motivated by 
Anderson's resonating valence bond (RVB) theory [2[ 
many investigators have proposed that the QSL and VBS 
phases can be explained within the correlated effective \- 
filled band scenario. For moderate Hubbard U, an AFM 




FIG. 1: (color online) Diagonal bond order B' and spin struc- 
ture factor S(tt, tt) versus U for t'=0A. (a) J' = (b) J' = J. 
Points are calculated at a spacing of AU—0.1 The J interac- 
tion strengthens the AFM phase. 



semiconductor-paramagnetic metal (PM) transition oc- 
curs with increasing frustration. D-wave SC mediated 
by fluctuations of the AFM ordering at the AFM-PM 
boundary has been also proposed based on mean-field 
and dynamic mean- field theories (DMFT) ■ 

Numerical calculations have, however, failed to find SC 
within the triangular lattice ^-filled band Hubbard model 
11- 3- Numerical studies have also failed to find a VBS 
phase in the model [141] . Although the ^ -filled Hubbard 
model on the anisotropic triangular lattice does not ap- 
pear to support SC, closely related models continue to 
be suggested as the appropriate theoretical model for de- 
scribing the SC transition in the CTS. It has been claimed 
that the simple Hubbard model does not include all the 
spin-spin interactions that play an important role in the 
CTS, and that additional spin exchange unrelated to the 
Hubbard U must be incorporated to correctly capture 
AFM fluctuation effects (Jo3]. This has led to theoret- 
ical works on the so-called Hubbard-Heisenberg model 
given below. The goal of this Letter is to critically exam- 
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ine whether the addition of a Heisenberg exchange term, 
assumed to be independent of the Hubbard onsite inter- 
action U, causes a superconducting phase to occur in the 
■i-filled band anisotropic triangular lattice model. A sec- 
ond objective is to see whether the combined effects of 
U, J, and frustration can mediate the VBS phase found 
in the EtMe 3 P[Pd(dmit) 2 ]2 [H- 

We consider the following Hamiltonian, 

H = -t -t'Yl Bi >i + U Y1 n *t™4 

<ij> {ij} i 

+ .1 ^ • j-y^s.-s,. (i) 

<ij> {ij} 

In Eq. [T] sites (ij) are nearest neighbors on a square 
lattice while sites {ij} are the next-nearest neighbors 
across a diagonal of each square plaquette. Sjj = 
So-( c l a c 3," + H.c), where c] a creates an electron of spin 
a on site i; — <y ia Ci a . Si is the (spin-i) spin operator 
for site i. All energies will be given in units of t. We 
limit our analysis here to the region of lattice anisotropy 
appropriate for the title materials, t' < 1 lil 18|. 

We consider two limiting cases: (i) J' = 0, when the 
added Heisenberg interactions do not frustrate Neel an- 
tiferromagnetism, and (ii) J' = J, which frustrates the 
AFM state. The J' — limit was studied by Gan 
et al using renormalized mean-field theory Q and the 
J' = J limit by Rau and Kee using a slave-rotor mean- 
field theory |10j. Powell and McKenzie studied variable 
J' [0 < J' I J < 2] using an RVB ansatz @. SC was 
found in some region of the parameter space by all of 
these authors. Several authors also found a spin-liquid 
phase 0, [T3] • Rau and Kee examined the isotropic lattice 
and in addition found a VBS phase [l(| that was claimed 
to explain the VBS in EtMe 3 P[Pd(dmit) 2 ]2 (H- 

Here we examine the ground state of Eq. [1] using ex- 
act diagonalization of a 4x4 lattice. Exact diagonaliza- 
tion was previously used to study the model in the limit 
J = 0, concluding that no SC or enhancement of the 
pairing correlations by U is present 11]. Despite the 



small lattice size (even smaller lattices are however used 
in cluster DMFT calculations Q) the validity of these 
4x4 results has been confirmed by recent Path Integral 
Renormalization Group (PIRG) [19( calculations on con- 
siderably larger lattices [12j that arrived at even stronger 
conclusions regarding the absence of SC. 

It will be useful to briefly recall the ground-state phase 
diagram in the J = J' = limit [TTI4l3. I20II2H . Known 
ground state phases include two AFM phases with Neel 
and 120° order, a PM phase and a sapl ess non-magnetic 
insulator (NMI) or QSL phase [HI 111 Q . The 

NMI phase is found between the AFM and PM phases 
for t' > 0.5 0, ED]. Near the isotropic lattice (*' » t) 
magnetically ordered states with Q ^ (vr,7r) are found 
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FIG. 2: (color online) Diagonal bond order B' versus U for 
t' = 0.8 with (a) J' = 0, and (b) J' = J. The derivative of 
B', dB'/dU with (c) J' = and (d) J' = J. B' is calculated 
for points separated by AU — 0.1. For J' — no NMI phase 
is found. When J' / 0, both PM-NMI and NMI-AFM phase 
boundaries shift to smaller U. The solid, dashed-dotted and 
dashed curves correspond to J=0.0, 0.3, and 0.5, respectively. 



We first consider the t' < 0.5 region of the phase di- 
agram where a direct PM-AFM transition is found for 
J = J' = 0. We have calculated the diagonal bond order 
B' = (Bij) for sites i and j connected by the t' bond 
and the spin structure factor 

s(Q) = ^E eiQ ' (rj_rfc) <Kt-^)Kt-^)>- (2) 

Fig. QJa) and (b) show B' and S(ir, n) versus U for t' 
0.4 and several different J (0, 0.3, 0.5). A sudden drop 
in B' at U = U c , simultaneous with a sudden increase in 
S(n, 7r), indicates transition to an insulating AFM state 



111 l21|. For both J' = and J' = J, J lowers the 



for large U [llj, [22J, [24[. However, such non-Neel AFM 
ordering is not found experimentally [l|. 



U c for transition to AFM, broadening the AFM region. 
As shown in Fig. Q] the largest broadening occurs when 
J' = 0. 

As t' increases the size of the discontinuity in B 1 and 
S(tt, 7t) decreases. This is due to the appearance of a NMI 
phase in between the PM and AFM phases for t' > 0.5 
[2l|. While the transition between PM and AFM 
phases for t' < 0.5 is discontinuous as U is varied, for 
larger t' the PM-NMI transition at U c i and the NMI- 
AFM transition at U c i are continuous [2l|. The transi- 
tions are now determined by examining the derivatives of 
order parameters with respect to U [llj. In Fig. [2] we plot 
B' versus U and its derivative dB/dU for t' = 0.8, calcu- 
lated using a centered-difference approximation with a U 
grid of AU = 0.1. We find that for J' = (Fig. Ufa) and 
(c)), the stability of the AFM phase is enhanced by U; 
U c i decreases with J. At the same time, the NMI phase 
is suppressed: the inflection point in B 1 in Fig. HJa) for 
J = and the corresponding minimum in dB'/dU in 
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FIG. 3: (color online) Long-range d x 2_ y 2 pair-pair correla- 
tion function P{r*) (see text) as a function of J for (a) t'=0A, 
J'=0 (b) £'=0.8, J'=0 (c) t'=0A, J'=J, and (d) t'=0.8, J'=J. 
P(r*) at U — J = is shown by the filled square on each 
plot. Independent of the value of J, long-range pair-pair cor- 
relations decrease monotonically as U increases. 



Fig. [H[c) disappear even for quite small values of J. For 
J' = J (Fig.[2jb) and (d)), the system progresses through 
two transitions, PM-NMI and NMI-AFM, as in the pure 
Hubbard model. As seen in Fig. [^b) and (d), there is a 
slight narrowing of the NMI phase. 

A necessary condition for SC is that the pair-pair cor- 
relation function P(r) for pairs of appropriate symmetry 
reaches a constant value as r — > oo. An additional re- 
quirement for SC mediated by interactions is that the 
pair-pair correlations are enhanced by the interaction 
[llj . In exact diagonalization studies of the J = model, 
P(r) for all symmetries was found to decrease monoton- 
ically with U from the [7 = limit [llj]. Large-lattice 
studies found further that as a function of distance r, the 
magnitude of P(r) decreased with distance faster than 
the U — solution [l2f . These two results indicated that 
SC is not present in the J — limit. Here we focus on 
any possible enhancement of P(r) due to J. 

We calculate the pair-pair correlation function as a 
function of distance, P(r) — (AjA,+f?), where pair- 
creation operators A- are defined as 



(3) 



The phases g(v) determine the symmetry of the super- 
conducting pairs. In our calculations we considered s, 
d x 2_ y 2, and d xy pairing symmetries [llj. Out of these 
pairing symmetries, we found that for if < t, P{r) is 
largest for d x i_ y i pairing symmetry. 

Although we have calculated P(r) for all r we show 
here our results for the largest possible r, r* . Fig. [3] 
shows the pair-pair correlation P(r*) as a function of J 
for d x 2_ y 2 pair symmetry. Representative values of U are 
chosen in each panel to correspond to the different regions 
of the phase diagram (PM, AFM, NMI). In Fig.^a) and 



(b) we take J' = and two different values of t', 0.4 
and 0.8. In both Figs. G^a) and (b]a direct PM-AFM 
transition exists. As with J = [11|, P(r*) decreases 
monotonically with U, with its value decreasing discon- 
tinuously at the PM-AFM transition. As Figs. Ela) and 
(b) show, P(r*) also decreases monotonically with J for 
J' = 0. The primary effect of J here is to expand the 
AFM region (Fig.[T]), which decreases the potential phase 
space available for SC. In Fig. EJc)-(d) we take J' = J. 
In Fig. (He), t' — 0.4 and again there are only PM and 
AFM phases. The behavior is very similar to the J' = 
case. Fig. |3][d) is for J' — J and t' = 0.8, with U values 
chosen to represent points in the PM, NMI, and AFM 
phases. In the PM region for large t' (if > 0.7) there is a 
small increase of P(r*) with J, but P(r*) remains consid- 
erably below its uncorrelated value. Pairing correlations 
continue to decrease monotonically with U. 

Summarizing our results for the d x 2_ v 2 pair-pair cor- 
relations, for t' < 0.7, J and J' decrease P(r). In the 
larger t' region, although there is a weak enhancement of 
P(r) compared to J = J' = 0, the pair correlations never 
exceed the value obtained for the uncorrelated limit. 

The compound EtMe3P[Pd(dmit)2]2 exhibits a spin 
gap below 25K [l5j]. The low temperature phase has 
been described as a VBS by the original investigators 
[l5| . Pressure induced VBS-SC transition, - analogous 
to the AFM-SC and QSL-SC transitions in the k-(ET) 2 X 
- occurs in EtMe3P[Pd(dmit)2]2- Any candidate model 
that is valid for CTS with dimcrized units should there- 
fore have VBS order in some region of the phase diagram. 
Further motivation to find VBS phases in models of inter- 
acting electrons comes from the extensive recent theoret- 
ical investigations of frustration-driven AFM- VBS quan- 
tum phase transitions within quantum spin models [25| . 
The VBS phase, if it at all appears within the present 
model, should appear in the highly frustrated region of 
the phase diagram. This is in agreement with the estima- 
tion of t'/t w0.9 in EtMe 3 P[Pd(dmit) 2 ]2, obtained from 
ab-initio calculation [18(. In the context of the present 
model the parameter region of interest is large U with 
if ~ t, in between the Neel and 120° AFM phases. Pre- 
vious PIRG numerical studies [3] as well as more recent 
work 26] on the Hubbard (J' = J = 0) model have, 
however, found the NMI rather than VBS in this region. 

As with SC, the VBS phase has been claimed within 
slave-rotor theory for the isotropic Hubbard-Heisenberg 
model (if = t, J' = J) for both J = and J > [l(|. 
The J interaction was found to strengthen the VBS order 
[Tol | ■ We therefore investigate the Hubbard-Heisenberg 
Hamiltonian in this highly frustrated region. The VBS 
order parameter is the bond-bond structure factor Sg(Q) 
for bonds along the x axis, defined as 
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FIG. 4: (color online) (a) Bond-Bond structure factor 
SBxilxiqy) as a function of q x , q y for t' = 0.9, U = 8, and 
J' = J = 0.3. (b) Sbx (w, q y ) for U = 8 and varying J. Lines 
are only guides to the eye. 
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FIG. 5: (color online) The phase diagram with J' = J for 
parameter regions relevant for the dimerized CTS. 



and either eliminates the NMI phase (for J' = 0) or re- 
duces its width in the phase diagram (for J' ^ 0). Impor- 
tantly, neither SC nor VBS phases are present in the |- 
filled band model for any realistic J. The phase diagram 
(Fig. [SJ remains nearly identical with t he p hase diagram 

The absence 



of the bare \ 
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filled Hubbard model 
of SC and VBS suggest that while the approximate meth- 
ods used previously can yield classical ordering such as 
AFM, spin-singlet formation, a quantum mechanical ef- 
fect, may be outside the capability of these techniques. 
The question that we started with - what is the mini- 
mal model that describes the CTS - then continues to 
be relevant. The only choice appears to be the inter- 
acting frustrated ^-filled band model [27| , which has the 
added advantage that it applies equally well to both the 
dimerized k-(ET)2X and Z[Pd(dmit)2]2 and the undimer- 
ized 0-(ET) 2 X which show CO-SC (as opposed to AFM- 
SC) transitions (recall that the charge carrier density per 
molecule is the same in these two classes of CTS). In 
recent work [3, [2t| we have shown that a frustration- 
driven AFM-spin singlet transition occurs within the 
dimerized interacting ^-filled band model, where the spin 
singlet state also exhibits CO. Structural analysis shows 
that the so-called VBS state in EtMe 3 [Pd(dmit) 2 ]2 has 



the same CO pattern [la, |30j as that in the coexist- 



CO-spin singlet found in the frustrated ^-filled band 
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. Experiments also find charge fluctuations in the 
QSL state of k-(ET) 2 Cu 2 (CN) 3 jUJ, suggesting again the 
applicability of the ^-filled band model. We are investi- 
gating whether frustration-induced SC occurs within the 
latter model. 

This work was supported by the Department of Energy 
grant DE-FG02-06ER46315. 



where (B) is the expectation value of the bond order. We 
consider just the J' = J case, as no NMI phase exists 
when J' = (see Fig. [5]). VBS order with columnar 
dimer pattern as claimed in 10] would correspond to a 
peak in S%(Q) at Q = (tt,0). Fig. gja) shows 5|(Q) 
for U = 8, t' = 0.9, and J' = J = 0.3. Throughout 
the NMI region (6 < U < 12) we find no changes in 
the Q-dependence of 5f>(Q). As shown in Fig. Q|a), all 
orderings with q x — 7r, corresponding to bond alternation 
along x, are suppressed compared to other values of q x . 
Fig. IDJb) shows the J dependence of S%(q x = Tv,q y ). As 
J increases, no peaks develop; rather the effect of J is 
simply a renormalization affecting all Sg(Q) at all Q 
equally. Since the VBS order is absent within the NMI 
phase for J = to begin with 1J] , this is strong evidence 
that the J interaction does not lead to a VBS phase. 

In summary, the inclusion of Heisenberg exchange in- 
teractions in the effective ^-filled band for for re^ET^X 
and Z[Pd(dmit)2]2 merely strengthen the AFM phase, 
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